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ABSTRACT 


This thesis contains the analysis of data of message 
traffic loads on the Naval Communications Fleet Multi- 
channel Broadcast System and results of a simulation model 
of the system under various channel alignments. Distri- 
butions of interarrival times, message lengths and requests 
for screens and eine are determined from the data. These 
distributions are used in a simulation model of the Broad- 
cast System. The model is used to compare the average 
delays caused by backlogs of messages when two channels 
devoted to a given ship-type are a) used in series with the 
second channel used to rebroadcast messages from the first 
channel after a one hour delay, and b) both channels are 
used in parallel to transmit first-run messages. The 
results of the simulation show that backlogs are reduced 
considerably by running the two channels in parallel at 


all times. 
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ΤΙ, FROBLEM 


This study concerns the message traffic loads on the 
Naval Communications Fleet Multi-Channel Broadcast System 
and the problem of determining the most efficient use of 
the 8 traffic channels in delivering teletype message 
traffic to the ships at sea. 

The Commander, Naval Communications Command message, 
R 0219442 SEP 71, officially commissioned this study. A 
copy of this message is shown in Appendix A for reference 
and to provide an example of a Naval Communications message. 

Appendix B provides background information on the 
organization of the Naval Communication System and a des- 


cription of the Fleet Multi-Channel Broadcast System. 


EX THE QUESTION AS POSED FOR SOLUTION 

The original question as it was posed for this study 
was: "When a significant backlog of any precedence message 
exists at a channel pair broadcast position, such that the 
primary channel delivering first run traffic is overloaded, 
when should the assigned secondary channel be activated for 
first run traffic and when that backlog ceases to be signi- 
ficant, when should the secondary channel be deactivated 
om first run traffic?". 

In order to discuss the problem further we must first 
consider how the primary and secondary channel pairs are 


aligned and used in different traffic load situations. 
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В, TYPES OF CHANNEL PAIR OPERATION 

Table XVII in Appendix B describes the fiset multi- 
channel broadcast alignment and the designated use for each 
of the 8 channels. This study considers only 6 of the 8 
— We note from Table XVII that these 6 channels 
form 3 pairs of channels which serve different types of 
ships copying the broadcast. They are: channels 1 and 2 
for the Destroyer force, channels 3 and 4 for the Service, 
Amphibious and Mine forces and channels 5 and 7 for the 
major warships. Channels 2, 4 and 7 are each designated 
for use either as a 1 hour delayed broadcast or as an over- 
load channel for the primary channels 1, 3 and 5 respectively. 
This allows these secondary channels to be employed in one 
of two different ways depending on the traffic load 
Situation. 

J. Situation 1 

When the traffic load for a channel pair ls such 

that the primary channel is not significantly backlogged the 
secondary channel is used for the 1 hour delayed rebroadcast 
of the primary channel's traffic. This alignment thereby 
ες the subscribers of the channel pair two chances to 
copy a message correctly and completely. Thus if a sub- 
scriber misses a message or a portion of a message on the 
primary channel the subscriber may then wait for one hour 
and attempt to copy the message on the rebroadcast channel. 
If the subscriber should miss the message both times he then 
must use the procedure described, in Appendix B, for obtain- 


ing the missed message. This configuration assumes that by 
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Broväding this redundancy in traffic delivery that the 
number of service requests for reruns will b reduced com- 
pared to having only one opportunity to copy the message. 
Further, it is noted that, for a subscriber to originate a 
screen request ts receive a reply and/or the desired 
rerun, the experienced normal time required for such action 
is in excess of one hour. Hence the 1 hour delayed rebroad- 
cast appears to provide better service proviced all sub- 
scribers use the system as this alignment intends. 

On the other hand there may be situations in which 
the urgency of the message is such that the subscriber may 
not want to wait and gamble on the 1 hour delay, only to 
possibly miss the message again, or perhaps equipment prob- 
lems may prevent a subscriber from being able to copy the 
rebroadcast. 

2. SAtüation 2 

This situation provides for the case when the 
channel pair traffic loads are heavy, and significant back- 
logs do build up such that a decision has to be made by the 
-NAVCOMMSTA to activate the secondary channel to transmit 
first run traffic in the same manner as the primary channel. 
In this case both channels of the pair are used to clear 
the backlog of traffic that has accumulated and to handle 
any anticipated increase in traffic arrivals which may 
cause any increased backlog. At the same time, however, 
this configuration precludes any facility for rebroadcasting 
the first run traffic as was described in situation 1 above. 


It is also assumed that this configuration will also prompt 
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more screen requests to be initiated by the subscribers 
since it intuitively seems reasonable that mcre messages 
may be missed since twice as much first run traffic is 
being delivered to the subscribers of the channel pair. 

In cases where a broadcast channel has no delayed 
rebroadcast facility use is made of what is called "heading 
recaps." Heading recaps are a periodic rebroadcast of the 
headings of messages previously transmitted on that partic- 
ular channel. Although the broadcast of the heading recaps 
will consume some of the valuable times needed for trans- 
mitting messages which are backlogged it is reasoned that 
this feature will have an effect of cutting down on the 
number of screen requests which were prompted because the 
Subscriber was trying to determine only if a missed message 
was of concern. 

Figure 1 diagramatically displays the above des- 
cribed situations, where channels 1 and 2 were the chosen 


channel pair for illustration. 


C. DISCUSSION OF THE PROBLEM 

Under conditions when no "significant" backlog of any 
message precedence exists the channel pair operation will 
operate as described in situation 1 above. However when a 
"significant" backlog starts to build up it may then be 
necessary to use situation 2 in order to make a more timely 
and efficient delivery of the traffic. 

The problem of this study is to determine what is meant 


by a "significant" backlog and thus determine the point at 
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Figure 1. Types of Channel Pair Operation. 
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which the channel pair iswswitched from situetion 1 to 2, or 
vice versa. At present this decision is mad by (i) observ- 
ing how the system has behaved in the past, (11) by observing 
the system in its present state, and (iii) by using experience 
to predict the ne state. 

The experienced communicator knows that the traffic 
arrival rate at the broadcast position follows a diurnal 
fluctuation. He has also observed when the channel pair 
operation was changed to the situation 2 ali;nment that the 
service desk after a time delay experienced a surge of screen 
requests since the rebroadcast channel facility was discon- 
tinued. This surge of screen requests which produce message 
reruns will most likely contribute further to the backlog 
and thus the backlog will get worse rather than better. 
Therefore when the backlog becomes significant the decision 
maker may decide to hold the system in situation l and hope 
that things get better, or he may make the decision to go to 
situation 2 and hope that the surge of screen requests. does 
not contribute to the backlog in a significant way. 

This method of decision making leaves something to be 
desired and this study was conducted in order to investigate 
the two alternate situations and to determine values of the 
decision variable for when switching should take place. 

When the NAVCOMMSTA activates the secondary channel for 
first run traffic, rather than as a rebroadcast channel, 
there is a 30 minute delay requirement between notification 


to the subscribers and the time the channel is activated. 
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This delay time used by NAVCOMMSTA San Francisco, and assumed 
common among other NAVCOMMSTA's, is needed f.r subscribers 

to prepare necessary equipment to copy the secondary channel 
all the time. Notification of the intent to activate the 
secondary channel is sent to the subscribers on the primary 


channel. 
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oe. Seen’ Oe ROACH TO THE PROBLEM AND ITS SOLUTION 


A. THE DECISION VARIABLE 

In order to answer the original question it is first 
necessary to determine when a significant backlog of messages 
exists. Let us say that a backlog is considered significant 
if the expected or average time a message of a certain 
precedence spends in the broadcast queue and transmission 
system exceeds the limits of the speed of service criteria 
shown in Table XVI in Appendix B. 

Actually the time standard referred to here is smaller 
than the speed of service criteria in the table, because the 
message has already spent some time in transmission before 
reaching the broadcast position. 

During the Peet cis of screen request messages in this 
study it was observed that these messages took an average 
ОҒ 2 hours in transmission from the originator to the ser- 
vice desk in the NAVCOMMSTA. If we observe that these 
screen requests are almost always of Immediate or Priority 
precedence and if their numbers were evenly divided then 
the average speed of service criteria would be approximately 
150 minutes. Thus 80 percent of the average transmission 
time was used in transmitting the screen request messages 
from the ship originators to the NAVCOMMSTA. If for example, 
these messages were not screen requests messages destined 


for the service desk, but rather for further relay on the 
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broadcast, we can see that on the average these messages 
could spend 30 minutes waiting at the broadc st position 
for transmission before exceeding the average speed of 
service criteria. We allow however that not all messages 
arriving at the Bere position are from ships at sea 
but that many are from shore establishments and hence the 
messages were received at the NAVCOMMSTA via the AUTODIN 
network which has a much faster transmission rate than the 
100 wpm ship-shore teletype circuits.  Never-heless it does 
not appear unreasonable to assume that an arbitrary but 
conservative revised speed of service criteria standard to 
be used at the broadcast could be 50 percent of the criteria 
shown in Table XVI in Appendix B. 

Reducing the times in Table XVI by 1/2 would result in 


times shown in Table I below. 


TABLE I 


BROADCAST SPEED OF SERVICE CRITERIA 


Precedence таты... Limits 
Flash (Z) As fast as Same 
humanly possible 
Immediate (O) 15 minutes 15-30 minutes 
Beiority (P) 1% hours 30 minutes-3 hours 
Routine (R) 3 hours 1%-12 hours 
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The decision variable of the problem can then be defined 
as the expected time a message of a certain »recedence 
spends in the broadcast queue and transmission system, i.e., 
average time in queue plus the average time in transmission. 
In addition, a measure of effectiveness can be defined to 
be the time difference between the expected time in system 
and the broadcast speed of service criteria. 

It would then be logical to conclude that if messages 
of a certain precedence spent more time on the average than 
is allowed by the broadcast speed of service criteria, the 
broadcast system is not making timely and efficient delivery 
of the traffic. Thus if this condition exists for any 
appreciable period it should be necessary to use both the 


primary and secondary channels to deliver the traffic. 


B. MODELING THE PROBLEM 

The type of model which best fits this system is a 
queueing model. The customers are the messages and the 
servers are the transmission channels. The system has four 
priorities or precedences of messages of which the highest 
precedence is Flash. The Flash message is allowed to pre- 
empt any lower precedence. When a flash message preempts 
another message the preempted message is later rerun from 
its beginning after the flash message transmission has been 
completed. A flash message is run three times to ensure 
receipt by the addressee. The preemptive method described 
is known as the Preempt Rerun method as opposed to the 


Preempt Resume method where the preempted messages will 
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continue transmission from the point at whick it was inter- 
rupted. 

The modeling of the feedback loop of this system is 
probably the most challenging because of its ne: time 
delays. The feedback loop referred to is that part of the 
system in which screen requests are made by the subscribers 
and from which the message reruns are generated. This feed- 
back loop has three time delays associated with it plus a 
message grouping effect and a change in precedence effect. 

The first time delay is the period of time between the 
actual time the message was missed and the time at which a 
screen request was sent out by the subscriber. This time 
delay was not analyzed due to the complexities of measure- 
ment. The second and third time delays are the time in 
transmission of the screen request to the NAVCOMMSTA and 
the time required to process the screen request at the 
NAVCOMMSTA broadcast service desk. 

The grouping effect is caused by the fact that screen 
requests normally require a number of messages to be 
screened and this in turn results in the rerun messages 
being delivered to the broadcast position in groups of 
messages. 

The change in precedence effect is something that has 
not been previously mentioned. It is observed that the 
normal precedence assigned to a screen request is either 
Immediate or Priority even though some of the messages 


asked for may be of Routine precedence. This is done because 
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the originator of the screen request senses a greater 
urgency in obtaining the missed messages because of the 
above described time delays. This upgrading of precedence 
then may require the NAVCOMMSTA to answer the screen request 
at a higher precedence than the precedence of any Message 


referenced in the request. 


C. REQUIRED: DATA 

The data required to analyze this problem falls into 
three areas. First, data to determine the type of message 
arrivals at the broadcast position. Second, the message 
length distribution must be known. Since the transmission 
rate of the servers is constant at 100 wpm, the message 
length distributions translate into the service time dis- 
tributions. Third, data concerning the screen requests is 
required. It is essential to investigate the time delay 
involved from the request for a rerun until it is received 
by the subscriber, as well as various other parameters which 


will be more fully described later. 


D. SUMMARY OF RESULTS 
In order to realistically model the problem, data for 
interarrival times of first run messages, message lengths, 
and requests for sceens and reruns was analyzed. The 
following conclusions were made. 
i. The first run message interarrival times were 
fitted to the exponential distribution and thus 


the arrivals follow a Poisson Process. (See 
Table II.) 
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ii. The message lengths were fitted to a hyperexpo- 
nential distribution with two prc ,ortion 
parameters. (See Table III.) 

iii. The interorigination times for screen requests 
fit the exponential distribution and thus were 
also concluded to arrive in a Poisson Process. 
(See Table X.) 

iv. The transmission delay distribution for screen 
requests was fitted to a right shifted expo- 
nential distribution. (See Table X.) 

у. The delay of a screen request at the service 
desk was concluded to follow right shifted 
hyperexponential distribution for high load 
periods and a right shifted exponential distri- 
bution for low volume traffic. (See Table X.) 

vi. The number of rerun messages resulting from a 
screen request was found to be exponentially 
distributed. (See Table X.) 

The relationship of the number of first run messages and 
the number of subscribers with the number of screens and 
reruns was analyzed by linear regression. The results indi- 
cate that these variables are independent of each other. 
(See Figures 12-19.) 

The GPSS simulation model of a broadcast channel pair 
was used to simulate the FASW channel pair in various con- 
figurations and traffic loads. The configurations simulated 
and compared by the average level of backlog were: i) channel 
2 closed to first run traffic and used as a rebroadcast of 
channel 1 traffic, 11) channels 1 and 2 both transmitting 
first run traffic while the rebroadcast facility was not 
used and 111) channel 2 alternately opened and closed to 
transmitting first run traffic whenever specified backlogs 
were exceeded. Тһе result of these simulations was that back- 


logs were reduced considerably by running the two channels 


in parallel transmitting first run traffic and reruns. 
23 





III. DATA COLLECTION AND ANALY.IS 


The primary source of data used in this thesis was 
NAVCOMMSTA San Francisco located at Stockton, California. 
The relative close location of Stockton to Monterey, Cali- 
fornia allowed frequent visits to the NAVCOMISTA both to 
talk to personnel at the station and to collect data. 

Although the fleet atmosphere of the Eas:ern Pacific 
differs from other parts of the world, it is felt that the 
type of data collected and that in general the analysis 
results are typical of any NAVCOMMSTA. 

All data which was analyzed to determine a goodness of 
fit to a theoretical probability distribution was analyzed 
by using the Fortran IV computer program described and shown 


in Appendix C. 


A. FIRST RUN MESSAGE INTERARRIVAL TIMES ANALYSIS 
P Source of Data 
The data used for the analysis represented 4 days of 
traffic received at the channel 1 broadcast position at 
NAVCOMMSTA San Francisco. The days represented were 5-8 
March 1971. The traffic load for these days was a low to 
medium volume of 142 to 401 messages per day. 
2. Analysis Results 
In order to show that the messages arrive in a 
Poisson Process it is necessary to show that the interarrival 
times are distributed according to the exponential prob- 
ШОНЕТ distribution. 
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The data was analyzed for the daily aggregated 
arrivals and by individual precedences. Tab.e II gives the 
results of the statistical goodness of fit tests. We see 
that the Kolmogorov-Smirnov tests indicated rejection of 
the hypothesis that the interarrival times are distributed 
exponentially while the Chi-square tests indicated acceptance 
for all samples except one. This apparent contradiction is 
explained by noting that the interarrival times were measured 
to the ες integer minute. Hence this caused the data 
to contain groups of data points of equal value, especially 
for times less than or equal to 3 minutes. As is noted in 
Appendix C this grouping effect causes the Kolmogorov- 
Smirnov (K-S) statistic to inflate and hence causes rejection. 

Figures 2 and 3 represent the graphic results for 
All precedences of 5 March, the data for which both statis- 
tical tests indicated rejection of the hypothesis. Figures 
4 and 5 represent the graphic results for All Precedences of 
6 March and are shown for comparison to the 5 March graphs 
since the chi-square test indicated acceptance of the 
hypothesis for the 6 March All Precedence data. 

From the statistical test results and from graphs of 
the data it was concluded that the interarrival times could 
be reasonably assumed to ke distributed according to the 
exponential distribution. Since this data represented only 
channel 1 arrivals, it is assumed that message arrivals at 


other channel pairs also arrive in a Poisson Process. 


` 
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B. MESSAGE LENGTH ANALYSIS 
1. Source of Data ste 

The data for message length analysis is from two 
sources. The first source contains two samples of messages 
and was taken at NAVCOMMSTA d Estes: These samples 
represented the traffic sent on 5 October 1971 on channel 1 
(FASW) and channel 5 (PNSC). "he other et of data was 
the Naval Electronics Laboratory Center. This data was 
obtained from the World Wide Traffic Demand Study data bank 
and it represents all messages sent on a fleet broadcast 
in various parts of the world on 8 October 1970. This 
sample therefore represented the Fleet Broadcast System as 
a whole and could be assumed not to be biased by any area 
effects. 

The FASW and FNSC samples were taken by measuring 
the messages in teletype tape form. This method of sampling 
provides for a high degree of accuracy in determining the 
length of a message. Teletype tape contains 100 characters 
per foot which equals 20 words. It should be noted that 
.the tape includes all teletype characters which are processed 
by the transmitter distributor when sending a message such 
as spaces, blanks, line feeds, shifts to upper and lower 
case print, etc., characters which are not easily identified 
on a printed copy. 

The message lengths of the World Wide sample on the 
other hand were determined bn silo page copies of the 
messages and using a regression equation to determine the 


number of words per inch of page copy. 
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2. Analysis of Message Lengths 


The general plan for this analysis wes to attempt 
to fit a probability distribution to the message length data 
by precedence and then to compare that fit to the world wide 
sample. 

The first fit attempted was the exponential distri- 
bution. The results of this showed that the message lengths 
were not distributed exponentially. 

At this point two observations of the above analysis 
were made. One was that a message has a minimum positive 
length, since each message has to have a heading which has 
a date-time-group, an originator and an addressee plus a 
text. It was estimated that such a minimum would be approx- 
imately 20 words long. This led us to fit the data with a 
right shifted exponential. However, the second observation 
that was made was that the graphs of the log of the recipro- 
cal of the tail distribution, or simply called the log of 
the tail, showed that the data in every case had a decreasing 
failure rate at least for message lengths up to approximately 
500 words. One particular distribution with a decreasing 
failure rate is a mixture of exponential distributions known 
as the hyperexponential distribution. 

These two observations led us to fit the data with 
a right shifted hyperexponential distribution with two 
parameters. 

The remaining question was whether the FASW and FNSC 


data which was taken from one day only was representative of 


32 





messages sent on other days and on broadcasts throughout 
the world. This is the primary reason that tne world wide 
sample was obtained. 

It was found that the world wide data behaved in 
the same manner as the FASW and FNSC data and that the right 
shifted hyperexponential made a very good fit. It was noted 
however that one of the major differences between the two 
sets of data was that the world wide data had very few data 
points which exceeded 1000 words while such cata points were 
very much in presence in the other two samples. Since the 
analyst personally took the FASW and FNSC samples it was 
decided that these large data points were not freak outliers, 
but were very much a part of the distribution. On the other 
hand it is noted that in reality, message lengths on the 
broadcast and in general do have an upper bound to their 
length which is estimated to be approximately 2000 words. 
This is an apparent contradiction between the empirical and 
theoretical distributions. It is resolved however by observ- 
ing that the probability of getting such a long message is 
indeed very small. 

3. Analysis Results 

Table III shows the results of the goodness of fit 
tests for message lengths. It should be pointed out that 
for all goodness of fit tests that not only were the statis- 
tical tests used in making a judgment on the fit but also 
interpretation of the log of the tail and the Cumulative 
Distribution Function (CDF) graphs. As is noted in Appendix 


C, the Chi-square statistic for large sample sizes tends 
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to inflate and careful attention must be made to the selec- 
tion of the number of intervals and their widths. In this 
analysis these drawbacks of the chi-square test were care- 
fully controlled. When the results of the K-S test and the 
chi-square test conflict we choose the K-S test results 
unless a graphical interpretation could explain the reason 
for the rejection result. 

We see in Table III that the K-S test results indi- 
cate that the E CREE should be rejected when all 
precedences are aggregated while only in the case of the 
FNSC data for all precedences does the chi-square test 
indicate rejection. Additionally the K-S test indicates 
rejection for Flash and Immediate precedences for the FNSC 
data. It is noted that the Flash message samples are 
included only to give representation of the parameters 
associated with this precedence since the respective samples 
were very small and give little meaning to any goodness of 
fit tests. | 

Figures 6 through 11 show the log of the tail dis- 
tribution and CDF graphs for the aggregated precedence for 
each of the three samples. 

In the cases where the Table statistical results 
indicate rejection for All Precedences of each of the three 
samples we see that the empirical curve in the log of the 
tail graph deviated from the theoretical curve in the message 
length range of greater than 600 words. However the CDF 


Eh Shows a reasonably good fit to the theoretical 
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distribution. Since the hyperexponential distribution is a 
mixture of two exponential populations one i terpretation is 
that the message length come from two exponential distribu- 
tions with the probability or proportion parameters p and 1-p 
as shown in Table III. It seems reasonable that all παρ 
lengths can be considered to be distributed according to the 


right shifted hyperexponential distribution. 


C. SCREEN REQUEST AND FEEDBACK LOOP ANALYSIS 
im Reasons 1or Analysis 

The analysis of the entire screen request and feed- 
back loop of this queueing system is a key part to the solu- 
tion of the problem. For if the system had no feedback then 
the solution of using both channels of a pair for first run 
traffic is obvious, since there would be no reason for 
rebroadcasting the first run traffic. 

The investigation of this area poses several new 
questions. First it is known that a certain percentage of 
messages are missed during the first run and for some sub- 
scribers these messages are successfully copied later on the 
 rebroadcast channel. In order to make some judgment of the 
effectiveness of the rebroadcast channel we can ask, "How 
many of the subscribers really do copy the secondary 
channel?" Next it is known that when messages are missed, 
that the screen requests usually contain a request for 
screening more than one message. This implies that the 
Subscribers do wait for a period of time before they ask 
for screens or they may miss several messages in succession 


or a mixture of both. This time delay was not analyzed 
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since it would be very difficult to collect data from the 
subscribers. 

When a screen request is originated there is a 
transmission delay between the origination time and the 
receipt by the NAVCOMMSTA. Also, when the screen requests 
are received by the NAVCOMMSTA there is a time delay for 
processing the request which is related to the number of 
messages that have to be screened and to the number of other 
screen requests waiting at the service desk. 

This entire process results in the feedback of 
messages to the broadcast. The distribution of rerun 
messages, times between them, and their relationships to the 
number of subscribers and first run messages must be 
analyzed. 

2. Source of Data 

The data for all of the screen requests and feedback 
loop analysis was obtained from NAVCOMMSTA San Francisco. 

3. Number of Subscribers that Copy Rebroadcast Channel 

A sample of 14 days of data was available for this 
analysis. During the period, 26 May-8 June 1971, subscribers 
were required to report to NAVCOMMSTA San Francisco periodi- 
cally as to whether they used the rebroadcast channel. Addi- 
tionally when a subscriber submitted a screen request he was 
to indicate whether he was using the rebroadcast channel. 
Table IV indicates the results in average numbers for each 


channel. 
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TABLE IV 


NUMBER OF SUBSCRIBERS USING 
THE REBROADCAST CHANNEL 


NR. SUBS 
REO SCRN 
NR. COPY το οσο 
NR. NR. NR. NR. NR.SUBS SECONDARY SECONDARY 


CHANNEL 1STRUN SUBS SCRNED RERUN REQSCRN CHANNEL CHANNEL 


FASW 293.1 14.6 26.2 : 6.5 3.9 EET 0.8 
(235485) (5.8%) 

FALD 243.4 19.5 24.8 5.4 Qut 2.6 0.4 
(13.532) (1.83) 

FNSC 165.9 11.8 39.2 5.1 4.8 3.2 2.0 


(27.25) (16.9%) 


Although this 1s a small sample and the area loca- 
tion for the sample does not represent the system as a 
whole it does however show that the secondary channel may 
not be used to its fullest benefit. If there is to be any 
conclusion drawn from these results it must be that a 
surprisingly small percent of the subscribers indicated that 
they used the rebroadcast channel. An explanation for this 
may be that the majority of the subscribers either do not 
desire to wait for an hour to attempt to copy the missed 
message and instead originate a screen request or they 
obtain the missed message from another ship. 

4. Analysis of the Causes of Reruns 

For this analysis two sets of data were available. 

One set covered the period of 26 May-8 June 1971 and the 


other the period of 25 August-17 September 1971. The 
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information in the sample contained daily stetistics on the 
number of first run messages, the number of . abscribers 
copying the channel, the number of screen requests processed, 
the number of messages required to be screened and the 
number of seien had to be rerun. This data repre- 
sented four different channels, FASW, FALD, FNSC and FRTT, 
the single channel fleet broadcast which is not part of the 
multi-channel system. 

The August-September data was broken into two parts. 
During the period 7-17 September a First Fleet Exercise 
бга ROPEVAL 3-71 was conducted and during this period the 
. Operations Evaluation Group (OEG) conducted an experiment 
of broadcast channel realignment on the FMUL broadcast [1,2]. 
During this period channel 4 was removed from use as a 
secondary channel to channel 3 and was converted into a 
dedicated exercise channel designated XRTT. ‘The XRTT 
channel carried only exercise traffic to the exercise 
participants and thereby lightened the load on the other 
channels during the exercise. This experiment had no direct 
relationship to this study except that the data collected 
during the pre-ROPEVAL period and the ROPEVAL period did 
provide for some comparisons. One such comparison was in 
the channel 3-4 configuration. Here the subscribers had a 
backup channel for rebroadcast for one period while during 
the exercise they did not and hence this would hopefully 
provide some insight into the effectiveness of the rebroad- 


cast channel. 
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The exercise time period also provided for moderately 
high traffic loads and a larger number of su»scribers. 
Additionally the exercise channel XRTT provided an oppor- 
tunity to look at a high volume of traffic being transmitted 
under simulated wartime conditions while this channel also 
had no backup channel for rebroadcast but used only heading 
recaps as did channel 3 during this period. 

Various statistics on screened and rerun messages, 
Called service statistics, are shown for the various 
channels during the August-September period in Tables V-IX 
and for the May-June period in Tables XVIII-XXI in Appendix 
D. These service statistics show no obvious relationships. 

In the service statistics Tables the Service Request 
Rate is abbreviated SVC REQ RATE, and is determined by the 


following equation: 


Nr. Messages Screened 


SVC REQ RATE = (Nr. First Run) (Nr. of Subscribers) 


This rate is used by some analysts of communications as a 
measure of feedback in the system which includes the rela- 
tionship of first run messages and subscribers [1]. Due to 
the great variation of this rate over the periods analyzed 
it was not clear whether there is in fact any relationship 
among these three measures. Therefore a simple linear 
regression analysis was used to determine if any such rela- 
tionship existed [3]. An attempt was made to determine if 


there was a linear relationship between the following: 
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Number of First Run Messages vs. Number of Screens 
Number of First Run Messages vs. Number ¿f Reruns 
Number of Subscribers vs. Number of Screens 
Number of Subscribers vs. Number of Reruns 
(Nr. of First Run) (Nr. of Subscribers) vs. Nr. of Screens 
(Nr. of First Run) (Nr. of Subscribers) vs. Nr. of Reruns 
These relationships were chosen, as in the case of 
the SVC REQ RATE, because intuitively one would suspect that 
an increase in the number of first run messaces broadcast in 
a day and/or the number of subscribers copying the channel 
would lead to an increase in missed messages and reruns. 
a. Results of the Statistical and Regression Analysis 
An important item in Tables V-1X, and those in 
Appendix D, is the percentage of first run traffic that 
result in reruns. This measure indicates the additional 
work that is placed on the system due to the feedback. We 
note that generally this percentage is between 2-4% even 
though the number of first run messages and the number of 
subscribers vary from channel to channel. Two exceptions 
are 8.2% for FNSC (Aug-Sep) and 18.2% for XRTT during 
ROPEVAL 3-71. Recall that ROPEVAL was an exercise simulating 
wartime conditions and hence this increase of reruns during 
this period is not surprising. Probably the explanation for 
this is that the subscribers sensed a greater urgency in 
obtaining missed messages plus they did not have a rebroad- 
cast channel associated with XRTT. This would tend to 


indicate that the rebroadcast channel may be an effective 
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feature of the system. Оп the other hand recall that during 
the ROPEVAL period channel 3 was operating w thout a rebroad- 
cast channel as a backup. Looking at the pre-ROPEVAL and 
ROPEVAL period for channel 3, in Table VI we see that the 
percent of reruns race 1.3% while the number of first run 
messages and the number of subscribers were also greater 
during the ROPEVAL period. This simple comparison appears 
to indicate that the presence of the rebroadcast channel had 
little effect on the number of reruns. An inportant appli- 
cation of the simulation model (described later) is to 
resolve these apparently conflicting observations. 

Figures 12 through 19 show representative graphs 
of the regression analysis. All graphs represent the 25 
August-17 September period. It is noted that the regression 
lines were not restricted to go through the origin even 
though it is obvious that if there were no messages trans- 
mitted there would be no feedback. Therefore the reader 
must be cautioned not to extrapolate the regression line 
beyond the range of the data shown. In interpreting these 
‘graphs one must consider the grouping or spread of the data 
points and the coefficient of the control variable. For the 
FASW graphs in Figures 12 and 13 we see that the regression 
line is nearly horizontal, i.e., it has a slope near 0. 
This indicates that the number of screens and reruns are 
independent of the number of first run messages times the 


number of subscribers. 
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The channel comparison mentioned above is shown 
in Figures 14 and 15. We note that although it appears that 
BEESe regression lines have significant slope the scale of 
the graph causes this and one should note again that the 
coefficients of the control variable are nearly zero and 
hence there appears to be no significant difference in the 
way channel 3 was Be 

Figures 16 through 19 provide results for the 
FNSC, FRTT and XRTT channels. We see that the XRTT regres- 
sion line has the largest slope of all the samples shown, 
however we must note that the sample size is nearly one half 
of the other samples. 

If any conclusions are to be drawn from this 
analysis it would be that the number of reruns are indepen- 
dent of the number of first run messages and the number of 
subscribers. However further analysis in this area appears 
necessary in order to fully accept this conclusion. 

Regardless of what conclusions can be made, 
these regression equations were used in the simulation model 
(described later) to mathematically simulate the effect of 
the rebroadcast channel on the system. 

5. Analysis of Time Delays in the Feedback Loop 
Since the model is in the form of a computer simula- 
tion it was required that each time delay of the feedback 
loop be analyzed and described for use in the model. 
The source for this data was the broadcast service 


logs from NAVCOMMSTA San Francisco. These logs contained 
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the following information: the ο πες of the screen 
request message, the time the screen request arrived at the 
service desk, the time the screen request processing was 
completed, the time an answer to the screen request was sent 
bsc on the broadcast, the number of screens each request 
contained and the number of reruns which resulted from each 
request. Logs from two different months were used. One 
represented a high volume traffic month which was September 
1971 and the other à low volume traffic month which was 
December 1971. 

Four time distributions were analyzed. These were, 
the interorigination times of the screen requests, the inter- 
arrival times of the screen requests at the service desk in 
the NAVCOMMSTA, the delay in transmission times for the 
Screen request messages from the originator to the NAVCOMMSTA 
and the delay in processing the screen request at the service 
desk. Additionally the distribution of the number of reruns 
per screen request was analyzed for use in the model. 


a.  Interorigination and Interarrival Distribution 
Analysis Results 


For each of these distributions the times between 
successive origination and arrival times were determined and 
then tested for goodness of fit to the exponential distribu- 
Exon by ae the computer analysis program. In the case of 
the origination time, the date-time-group of the screen 
request was considered as the origination time of the request. 
In both cases the hypothesis that these times were dis- 


tributed exponentially was accepted. The results are shown 
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in Table X and Figure 20. This result implies that screen 
requests arrive at the service desk accordir to the Poisson 
Process in the same manner as first run traffic at the 
broadcast position. Only the interorigination distribution 
is used in the se 

b. Transmission Delay Distribution Analysis Results 

The delay in transmission times were considered 

to be the time difference between the origination time and 
the time of arrival at the broadcast service desk. This 
data was tested against a right shifted exponential for 
goodness of fit. The right shift was used because of the 
observation that a message being transmitted from a ship to 
the NAVCOMMSTA has at least some minimum time delay. Table X 
shows the results for both the regular and right shifted 
exponential fits. The hypothesis that the data fit the 
right shifted exponential was accepted. Figure 21 shows the 
CDF graph for the September data sample. It is interesting 
to note that the mean transmission delay for the two sets of 
data were nearly the same. The sample mean for the December 
. data was 119.6 minutes and for September it was 120.8 
minutes. 

С. Service Desk Delay Distribution Analysis Results 

The service desk delay times were considered to 

be the time difference between the time of receipt or 
arrival of the screen request message at the service desk 
and the time the reply to the screen request was completed. 


It should be noted that the times measured here included 
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both the time the screen request spent in queue at the desk 
and the time for processing at the desk. Th.s was purposely 
done in order that this distribution could be used in the 
model to implicitly simulate the delay resulting from other 
screen requests for other broadcasts and/or channels being 
handled at this same desk. 

This data was tested against a right shifted 
distribution. The right shift was again used by the obser- 
vation that each request however small must take some time 
to be processed. Table X and Figures 22 and 23 show the 
results of these tests. In the case of the December data 
the right shifted exponential was accepted as the underlying 
distribution while for the September data, a right shifted 
hyperexponential was used for the best fit. 


d. Reruns per Screen Request Distribution Analysis 
Results 


The data for this analysis represented the 
number of 'reruns per screen request for all channels pre- 
viously listed for the August-September period. The data 
was fitted with a regular exponential distribution and the 
hypothesis was accepted that the reruns per screen request 
data is distributed exponentially. The results of the test 


are shown in Table X. 
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Figure 20. CDFs for Screen Request Interorigination 


Times (September). 
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Figure 21. CDFs for Transmission Delays (September). 
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Figure 23. CDFs for SErvice Desk System Delays (September). 
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IV. THE MODEL 


A. GENERAL DESCRIPTION 

The model is in the form of a General Purpose Simulation 
System/360 (GPSS) computer program. GPSS is a simulation 
system designed for simulation of queueing systems and hence 
adapts very well to simulating the fleet broadcast system. 


References 4 and 5 fully describe GPSS. 


B. MODEL CAPABILITIES 

The model generates transactions which are the messages 
in the system and moves these messages through the model for 
queueing and processing in a very similar manner as the real 
world broadcast system. Since each of the channel pairs of 
the broadcast operate EEE, of the others, the model 
was designed to simulate only one channel pair of the fleet 
multi-channel broadcast system. 

The model nes essentially 3 basic components or routines 
which interact with each other. They are (i) the first run 
traffic flow for Immediate, Priority and Routine precedences, 
(ii) the Flash message routine with preemption, and (iii) 
the screen request and feedback loop routine. For the 
description of the model we will use "channel 1" to refer to 
the primary channel and "channel 2" to refer to the secon- 
dary channel. 

The model is essentially very simple in its flow charac- 


teristics and has much flexibility in its use. This 
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flexibility is described by listing the following features 
which can be considered as measures of controlling the model. 
The user may change or specify the following control vari- 


ables, distributions or methods listed. 


i. The daily number of first run messages to be 
generated for the simulation period. 

ii. The percent of each precedence of the first run 
messages generated. 

iii. The individual message length distributions for 
each precedence. 

iv. The mean message length for each precedence. 

v. The diurnal message arrival function. 

vl. The number of subscribers. 

vii. The backlog levels by precedence to be allowed 
before opening and closing channel 2 for trans- 
Mistting first run traffic. 

viii. Allow Flash messages to preempt other messages 
using either the Preempt Rerun method or the 
Preempt Resume method. 

ix. The mean transmission delay for screen requests 
or the distribution of this delay. 

x. The percent of screen request to be processed 
by the Broadcast Control Station and the percent 
to be processed at other NAVCOMMSTAs. 

X1. The service desk processing delay distribution. 

X11. The mean service desk delay for the BCS and/or 
the other NAVCOMMSTAs. 

xiii. The mean number of reruns per screen request or 
the distribution. 

xiv. The mean transmission delay for reruns to reach 
the BCS after being processed by the other 
NAVCOMMSTAs. 

XV. The control of whether rerun messages processed 


by the BCS are allowed to go to the head of the 
line of first run messages waiting in queue. 
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Ci DESCRIPTION Or THE MESSAGE FLOW IN THE MODEL 

The flow charts shown in Figures 24 through 27 provide a 
simplified picture of the message flow through the model. И 
First consider the traffic flow for Immediate, Priority and 
Routine precedences. These messages plus Flash are generated 
according to an exponential distribution with a mean inter- 
arrival rate specified and a specified diurnal function. 

After a message is generated its precedence is determined by 
a percentage function specified. After the precedence is 
assigned the model routes the messages by precedence to be 
assigned their appropriate message lengths from functions 
specified and then puts the messages in queue by precedence. 
At this point in the model the Flash messages are handled in 
a slightly different manner than the other messages. This 
is described below. The precedences or priority numbers 
assigned by the model, depending on the precedence assignment 
function specified, are as follows (listed from lowest to 
highest with designations in parentheses) 

1-Routine First Run (R) 

2-Routine Rerun (RRR) 

3-Priority First Run (P) 

4-Priority Rerun (PRR) 

5-Immediate First Run (0) 

6-Immediate Rerun (ORR) 

7-Flash (First Run and Rerun) (z). 

The model chooses the highest precedence and in order of 
arrival within precedences when advancing transactions through 
the model. Thus messages waiting in queue are chosen by 
highest precedence and on a first in first out basis for 


routing to the transmission channels. We note that when 


priority numbers 2, 4 and 6 are assigned to reruns processed 
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by the BCS, these messages will be transmitted before their 
first run counterparts. Hence this simulates allowing the 
reruns to be put at the head of the line. 

When a message is routed for transmission the model 
checks to determine the number of backlog by precedence and 
compares it to a specified value of allowed backlog for a 
given precedence. If the backlog for any precedence exceeds 
the allowed backlog the model will start a 30 minute timer 
(or a time specified) which represents the notification to 
the subscribers of the intent to activate channel 2 for 
first run traffic. At the end of the 30 minute period 
channel 2 is opened for first runtraffic. Whether or not 
the timer is activated the message proceeds to channel 1 
for transmission. ТЕГЕН the flow chart best illustrates 
this test as occurring sequentially between the time the 
message departs the queue and before transmission, this is 
not the way it occurs in the model. Actually all the above 
described assignments and tests are made at the instant in 
simulation time when the message is generated. That is, 
when the message is generated even though it is part of the 
queue, the message proceeds through the model until it is 
refused entry to some GPSS block (program statement). Thus 
the messages proceed through the GPSS model to the BUFFER 
block and the main GATE 1 before being stopped. For statis- 
tical purposes the messages do not depart from their respec- 
tive queues until they pass through the DEPART block and 
this occurs only when a channel is available and the message 


is allowed to SEIZE a channel in order to be transmitted. 
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The message is held in the transmission channel (ADVANCE 
block) for a period ος time equal to its assigned message 
length. When the transmission is complete the message 
releases the channel and is terminated and removed from the 
system. The main GATE 1 closes when a message SEIZEs a 
channel and opens when the message RELEASEs a channel. 

If channel 2 has been activated for first run traffic, 
then the model processes the messages through channel 2 in 
the same manner as in the primary channel. Messages waiting 
for transmission are selected in the precedence order 
described above and each time a channel is finished with a 
transmission, or is available, the model will move the 
available message to the appropriate channel. When a mes- 
sage is processed by channel 2, the backlog of the Priority 
and Routine precedences are checked against a specified value 
to determine if channel 2 should be closed. If channel 2 
can be closed this will activate a timer which runs for the 
duration of the message length of the message in channel 2 
and then closes channel 2 to any further first run traffic. 

It is noted with emphasis that the secondary channel in 
the model does not process messages when it is considered as 
being used for rebroadcast of the primary channel traffic. 
In the model this feature is handled through a mathematical 
relationship which will be more fully described below. 

The next N of the model to be described is Ehe Flash 
message routine. This is described separately to point out 


that although the Flash messages are generated as a part of 
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the overall first run traffic described above, these messages 
follow a different route through the model because they are 
allowed to preempt other lower precedence messages. When a 
Flash message arrives in the system it immediately preempts 
any message in channel 1 and begins transmission. Just 
before transmission and in the simulation time instant of 
generation the Flash message is assigned its message length 
along with a 25 word header and then the entire length is 
multiplied by three. This part of the simulation represents 
two things. First before a Flash message starts transmission 
on the broadcast channel a heading is sent which reads, 

"ZUJ FOR FLASH---------ZUJ FOR FLASH" etc. , and means Stand 
By for a Flash Message. This header also causes attention 
bells to ring on the subscriber's teletype thereby drawing 
attention to the fact that a Flash message is about to be 
transmitted. The second item is that a Flash message is 
always transmitted three times in succession to ensure re- 
ceipt by the subscribers. 

The Flash routine also provides for the situation when 
more than one Flash message may arrive. In this situation 
if channel 2 is being used for first run traffic and chan- 
nel 1 is already processing a Flash message then the second 
Flash will preempt any message in channel 2. 

When a message is preempted by a Flash message, the 
preempted δες is handled by the model in such a manner 
that when the Flash finishes transmission, the preempted 


message will be run again from its beginning, if the model 
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is being operated in the Preempt Rerun mode. If channel 2 
is open for first run traffic and a message is preempted in 
channel 1 and if channel 2 becomes available before the 
Flash has finished transmission in channel 1 the preempted 
message will be routed and transmitted on channel 2. 

The third part of the model is the Screen Request and 
Feedback Loop of the broadcast system. Screen requests are 
generated according to an exponential interorigination dis- 
tribution with- specified mean value. When the screen 
requests are generated they are assigned a transmission time 
from a right shifted exponential distribution with a given 
mean time. This holds the request in an ADVANCE block for 
this period of time before allowing it to proceed through 
the model. It is noted that more than one request can be 
in this ADVANCE block at one time thus simulating any number 
of screen requests in transmission to the NAVCOMMSTA. The 
model then randomly determines according to a percent speci- 
fied whether the screen request will be handled by the BCS 
or by another NAVCOMMSTA. Next the request is assigned a 
service desk system processing time and is held in an ADVANCE 
block to simulate the time a screen request spends in the 
service desk queue and on the service desk being processed. 
This time delay distribution with its specified mean repre- 
sents the competition a screen request for the channel pair 
being simulated has with other screen requests for other 
broadcast channel pairs or other broadcasts. When the 


request has finished processing and is ready to proceed, 
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Figure 27. Screen Request and Feedback Loop. 





it is then assigned the number of reruns that have resulted 
from this screen request from an exponential distribution. 
If the number of reruns is equal — the screen request 
transaction is terminated and removed from the system. If 
the number of reruns is greater than or equal to one, then 
each rerun becomes a message itself and is assigned a prece- 
dence and an appropriate message length and sent to its 
appropriate queue to await transmission. The option of 
allowing reruns to go to the head of the line was described 
above. If a screen request was processed by a NAVCOMMSTA 
other than the BCS, those reruns are assigned a transmission 
time from the transmission delay distribution with a speci- 
fied mean and held for that period in an ADVANCE block before 
joining their respective queues. This time.simulates the 
time it takes for the reruns to be transmitted via AUTODIN 
to the BCS. It is noted that in this case the reruns cannot 
allowed to go to the head of the line since these rerun 
messages arrive at the broadcast position in the same manner 
as first run messages and therefore are handled in the same 


Manner. 


D. THE MODEL AND THE REBROADCAST CHANNEL SIMULATION 

As was mentioned above, this model does not rebroadcast 
first run messages as though they were being sent on al 
hour delay rebroadcast channel. The secondary channel in 
the model is used solely for processing first run messages 
and when it is not doing so it is not used. The effect of 


the use or non-availability of the rebroadcast channel to 
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the subscribers is simulated through the control of the mean 
interorigination times of the screens requests. This is 
done by using a representative regression equation of (Nr. 
First Run) (Nr. Subscribers) vs. Nr. Screens for the channel 
pair being simulated. Thus by specifying the number of 
first run messages and the number of subscribers for the 
simulation period, the model determines the number of 
screens that will be handled. Then by further specifying 
the average number of screens per request the model deter- 
mines the number of screen requests to be generated and with 
this information the model calculates the mean interorigina- 
tion time by dividing the number of requests into the 
simulation time period. In the real world system it is 
assumed that when the керк ct channel is being used, 
that the number of screens requests is less than when the 
secondary channel is not being used for rebroadcast of 

first run traffic. Then if this assumption is in fact true 
the analyst need only decrease the mean interorigination 
time for the screen requests when the secondary channel is 
open for first run traffic. This therefore simulates a 
higher number of screen requests being generated and handled 
at the service desk and thus more reruns are generated. 

This is done in the model by one of two ways. The regres- 
Sion equation used can be changed when channel 2 1S open or 
Simply the number of screens per request can be changed in 
Order to effect the increase of screen requests. Note that 


this change in the mean interorigination time is made at 
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the time that the secondary channel is opened or closed for 
Erst run traffic. Thus we can observe how the increase in 
the number of reruns affect the system and since the time 
delays are built in the feedback loop we can also observe 
any surge in reruns which may occur after the secondary 


channel was opened for first run traffic. 


E. VERIFICATION OF THE MODEL 

As the model was built, each routine added was run 
separately to verify that the logic design was correct and 
the transaction flow was being performed correctly. The 
means and variances of interarrival times and message 
lengths, etc., generated by the model using random number 
generators and the specified functions, were checked for 
accuracy. Appropriate modifiers (explained in model docu- 
— were chosen to ensure that the model generated 
the function values with a reasonable accuracy for the mean 
values desired. 

One important verification was made by running the model 
for the M/M/1 and M/M/2 queueing systems. That is, the 
model was run without the feedback loop and the Flash pre- 
emption routine and with exponential (M) interarrival and 
service times for one server (only channel 1 open) and two 
servers (both channels open) [6,7]. The results of these 
model runs were compared to analytical calculations for the 
average server utilization and the average queue length and 
waiting times. These comparisons showed that the model was 


very accurate. 
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The model was also verified in its final form by using 
data obtained from OEG which was gathered during ROPEVAL 
3-71 [1]. This data represented the actual hourly number 
of arrivals and backlogs observed during the period of 
13-16 September on channel 1 (FASW). For the comparison 
runs with the model, parameters of the model were set at 
values that were observed during this time period. The 
parameters which were not varied for the comparison runs 
and for later described runs are shown in Table XI below. 

The observed traffic loads, number of subscribers number 


of screen requests and number of reruns are shown in Table XII. 


TABLE XI 


MODEL PARAMETERS HELD CONSTANT 


Flash Immediate Priority Routine 


Msg Lgth 89; 2 | 2453 25753 182.6 
Mean 

$ of First 1529 2245 3557 40.3 
Run 

$ of Reruns 1725 2s. 5 35% / 40.3 

Allowed 

Backlog for 
Opening 3 10 50 110 
CHAN2 
Closing _ ~ 1 2 
СНАМ2 
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TABLE XII 


FASW 13-16 SEP OBSERVED TRAFFIC LOADS 


Nr. Nr. 

Nr. Nr. Screen Msgs 

Day First Run Subscribers Requests Rerun 
13 SEP 422 38 6 6 
l4 SEP 668 40 8 4 
15 SEP 598 40 13 10 
16 SEP 480 33 6 36 
Totals 2168 33 56 


To simulate the daily and hourly fluctuations of traffic 
arrivals during the period 13-16 September the model used 
the number of first run messages shown in Table XII above 
and the "diurnal" function shown in Figure 28 below. The 
specified dirunal function allows the model to determine 
the proportion of a day's first run traffic that will arrive 
or be generated each hour of the simulation run. Another 
function in the model specifies the number of subscribers 
copying the channel pair each day, the numbers used are 
shown in Table XII above. The percentage of first run 
messages and reruns that were assigned each precedence are 
shown in Table ΧΙ above. These percentages were the observed 
aggregate percentages for the 13-16 September period being 


Simulated. 


86 





935 91 
р0 06 91 οἱ 


-SsTeaTıay 9ρεςς59θ MSV4 "AUS 9T-ET 201 чотаоЧод теыхата 


(SWIL NINZ ) Ava YO ЗМІ 


435 ΘΙ 
8 t vo OC 9| el 


8 


vb bd OS 


91 


335» δι 
οἱ 


8 


t 


82 әхпбтя 


djs εἰ 
ve OO 9I 


οἱ 


O 


¿0 


РО 


99 


80 


SIVAIJAV NNY 158414 
SAVA JO NOI.LJOdOHd ATYNOH 


87 





In the simulation the Broadcast Control Station (BCS) 
handled all screen requests. The parameters used for the 
feedback loop were: 

i. Mean Number of Screens per Request = 8 
ii. Mean Number of Reruns per Request = 1 


iii. Mean Screen Request Transmission 
Delay = 120 mins. 


iv. Mean Service Desk System Delay = 75 min. 
To determine the number of screens (messages to be screened) 
that would be handled each day the regression equation, 
Y = 44.36 + .00065 (X), was used, where X is the number of 
messages transmitted the previous day times the number of 
subscribers. This equation is taken from the above described 
regression analysis for FASW (channel 1) during the 
25 August-17 September period. Recall that during this 
period channel 2 was used for rebroadcast of channel 1 
traffic. The model determines the mean interorigination 
time for screen requests by first calculating the number of 
screens and then dividing that number by the specified mean 
number of screens per request. This gives the number of 
screen requests. By dividing the number of screen requests 
into 1440 minutes the result is the daily mean interorigina- 
tion time for screen requests. 

To compare the model results to those actually observed 
during the period 13-16 September the model was run in eo 
different forms. First the model was run with empirical | 
message length functions and with empirical functions for 


the service desk system delays and the number of reruns per 
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screen. Secondly, the model was run with theoretical func- 
tions for all distributions except Flash message lengths. 

The Flash message length distribution was left as an 
empirical distribution since the size of sample tested did 
not allow us to conclude the exact nature of its distribution. 

The results of these comparison runs are shown in Table 
XIII and Figure 29 below. By comparing the average queue 
lengths and the average wait in queue we see that the model 
is conservative. We also note that the model with theoreti- 
cal functions is more conservative than the run with empirical 
functions. In comparing the backlog curves in Figure 29 we 
note that the model tracks well. That is, the model shows 
peak backlogs build ups and declines at nearly the same times 
as were observed from the actual data. 

It should be pointed out that we cannot say that this 
comparison validates the model since each simulation run and 
the actual observations for this period are observations of 
a random process. Even if the actual data could be observed 
again in the very same circumstances as the period of 13-16 
September the results would not be the same, nor would 
Simulation model runs with different seeds for the random 


number generators be the same. 


F. MODEL SIMULATION RUNS AND RESULTS 


1. Description of the Simulation Runs 


The model with theoretical functions was used for 


this experiment since it was determined above that it was 
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more conservative than the model with empirical functions. 
The channel pair —" was simulated under three different 
situations, i) channel 2 closed to first run traffic but 
considered as being used for rebroadcast of channel 1 
traffic, ii) channels 1 and 2 both transmitting first run 
traffic and the rebroadcast facility is not used, and iii) 
channel 2 is opened and closed during the run by the timers 
when specified backlogs were exceeded. Additionally the 
above situations were run where reruns were allowed to go to 
the head of the queue as opposed to going to the back of 
the line. 

Eight simulation runs were made with these config- 
urations. Four simulation runs were made representing the 
period of 13-16 September and the other four runs represented 
a hypothetical four days in which the first run traffic loads 
were increased by 200 messages per day over those observed 


in Table XII above. The runs made are listed as follows: 


Run Iw 19-16 SEP. Channel 2eChosedefor First Run 

Run 2. 13-16 SEP. Channel 2 Open Whole Period 

Run 3. 13-16 SEP. Channel 2 Opened and Closed by Timers 

Run 4. Same as Run 2 with Reruns going to the back of 
the line 

Run 5. High Loads, Channel 2 Closed for First Run 

Run 6. High Loads, Channel 2 Open Whole Period 

Run 7. High Loads, Channel 2 Opened and Closed by Timers 

Run 8. Same as Run 7 with Reruns going to the Head of 


the Line. 

In runs 1-3, reruns were allowed to go to the head 
of the queue and for runs 5-7 reruns were put at the end of 
the queue. All of the runs used the dirunal function shown 
in Figure 28 to determine the hourly first run arrivals. In 


runs 2-3 and 6-7, when channel 2 was open, the number of 
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screens were calculated by using a different regression 
equation and mean number of screens per request such that 
the model would simulate a higher number reruns. For runs 
1-4 the mean service desk system delay set at 75 minutes 
while for runs 5-8 it was set at 158 minutes (see Table X). 

The allowed backlog levels shown in Table XI for 
opening channel 2 in simulation runs 4 and 7 were chosen 
based on the Broadcast Speed of Service Criteria shown in 
Table I. By using a familiar queueing theory equation, 
L = AW, where L is the average queue length, W is the average 
wait in queue and A is message arrival rate, one can deter- 
mine the approximate backlog levels allowed before the Speed 
of Service Criteria are exceeded [6,7]. We note that this 
equation relates averages over a period of time. Therefore 
the maximum allowed routine queue length, which is usually 
the largest queue, was determined by using this equation and 
the arrival rate of routine messages observed during the 
13-16 September period. The other allowed backlogs by 
precedence were chosen somewhat arbitrarily but such that 
the allowed backlog would not exceed on the average, the 
Broadcast Speed of Service Criteria. 

2. The Results of the Simulation Runs 

The results of these eight simulation runs are shown 
in Tables XIV and XV and Figures 30-32. In comparing the 
average queue lengths, average wait in queue and the maximum 
backlog observed plus the backlog curves the result of the 


simulation is obvious. The system configuration where both 
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channels transmit first run traffic and the rebroadcast 
facility is not used is or to any other configuration. 
Recall that the decision variable Bur this problem is the 
average wait in queue of any precedence message such that 
this time does not exceed the Broadcast Speed of Service 
Criteria indicated in Table I. We also note that the primary 
goal of the fleet broadcast is to make timely and efficient 
delivery of traffic to all subscribers. The model shows 

that there is no question as to the timeliness of traffic 
delivery when both channels are running first run traffic. 

So now we must question whether this configuration is also 
the most efficient method. This can be answered by looking 
at the feedback loop and the message reruns. oe the 
results for the high load runs 5 and 6 where channel 2 was 
closed for first run messages in one run and open for first 
run in the ewe солан For these two runs the mean time delay 
in transmission of a screen request was 120 minutes and the 
mean service desk system delay was 158 minutes thus on the 
average it took 278 minutes (4.6 hours) for a screen request 
to reach the broadcast position for transmission since the 
time it was originated. Now if we chose the routine messages 
for comparison, we note that when channel 2 was closed to 
first run traffic the average time it took for a routine to 
be rerun is 2433 minutes (40.6 hours) as compared to 326 
Minutes (5.4 hours) when channel 2 was transmitting first run 
traffic. We also note that 15.5% of the total first run 


traffic was rerun in run 6 as compared to only 4.4% in run 5. 
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To fully explore this area we further look at the 
comparison of the runs where reruns were allowed to go to 
the head of the line. We see that although the reruns did 
not have to wait as long in queue as they did when they went 
to the back of the line, we also see that the first run 
traffic had a longer wait in queue on the average than when 
reruns were put at the back of the line. This leads us to 
ask the question, "Which type of message is most important?" 
It would appear that the first run message is most important 
because it has not been transmitted to any of the subscribers 
copying the broadcast and hence it is really undisclosed 
information while a rerun message has most likely been 
received by the majority of addressees, if more than one, 
and hence it is old information to those subscribers except 


fem the subscriber who missed the message. 


100 





V. CONCLUSIONS 


The conclusions which can be made are that the broadcast 
channel pair configuration where both channels transmit first 
run traffic is the superior configuration for timely delivery 
of traffic. The comparison of time delays for a subscriber 
to receive a rerun show that this configuration is not only 
timely, it is efficient as well. Therefore, the answer to 
the problem question is that the system should always be run 
with the parallel channels transmitting first run traffic 
and that the rebroadcast of first run traffic on the secondary 


channel should be discontinued. 
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APPENDIX A 


ORIGIN OF THE STUDY 


. The problem on which this study is based was proposed as 
a thesis topic to this student in April 1971 by Captain F. M. 
Snyder, USN, then the Assistant Commander for Operations and 
Readiness to the Commander, Naval Communications Command 
(COMNAVCOMM). The study was officially commissioned as a 
thesis study by COMNAVCOMM message R 0219442 SEP 71. A copy 
of this message is shown below for reference and as an 


example of a Naval Communications teletype message. 


С2ССААЙ50 

RTTUZYUW RHELSAAQ965 2452994-UUU--RUWJAGA 

ZNR UUUUU 

R 8219442 SEP 71 

FM COMNAVCOMM 

TO RUWJAGA/NAVAL POSTGRADUATE SCHOOL 

INFO RUWNSAA/NAVCOMMSTA SFRAN 

BT 

UNCLAS //Ngg9ggg// 

FOR LCDR. HARLAN D. OELMANN 

THESIS STUDY 

1. YOUR THESIS STUDY OUTLINE HAS BEEN REVIEWED AND ITS 
PURSUIT IS STRONGLY ENCOURAGED. COMNAVCOMM REQUESTS 
ESETES OF THE COMPLETED STUDY AND ANY INTERIM REPORTS 
THAT ARE DEVELOPED. 

2. DR. HARDY (OEG) WILL VISIT NAVCOMMSTA SANFRAN DURING 
EARLY SEPT., AND IT IS SUGGESTED THAT YOU MEET OR SPEAK 
WITH HIM AT THAT TIME IN ORDER TO BECOME APPRISED OF 
OEG EFFORTS IN THIS FIELD. 

3. SHOULD COMNAVCOMM ASSISTANCE APPEAR NECESSARY OR 
BESTRABTE IN THE DEVELOPMENT OF THE STUDY, REQUEST YOU 
MAKE YOUR NEEDS KNOWN VIA NAVCOMMSTA SANFRAN. 

BT 
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APPENDIX B 


PROBLEM BACKGROUND 


I. A BRIEF EXPLANATION OF THE NAVAL COMMUNICATIONS SYSTEM 


A. ORGANIZATION OF THE NAVAL COMMUNICATIONS SYSTEM 

The Naval Communications System is composed of 27 Naval 
Communications Stations (NAVCOMMSTAs), 4 Naval Communications 
Units (NAVCOMMUs) and several Naval Radio Stations 
(NAVRADSTAs) [8]. These stations are located throughout the 
world and are organized into several Naval Communications 
Areas (NAVCOMMAREAS). The primary purpose of these stations 
15 to relay messages to ana from ships of the fleet and 
shore establishments within their area of responsibility. 
The stations within a NAVCOMMAREA are organized under one 
NAVCOMMSTA which is designated as the Naval Communications 
Area Master Station (NAVCAMS) while the other stations are 
designated as Naval Communications Area Local Stations 
(NAVCALS). Although the communication means which each of 
these stations use take several forms, such as teletype, 
voice and CW, etc., this study is concerned only with tele- 
type communications and the fleet broadcast system. There- 


fore discussion will consider only those areas. 


В. TELETYPE COMMUNICATIONS 
It is assumed that to most readers, teletype communica- 


tions is a familiar means of communication which does not 
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require detailed explanation, except to note that when 
messages are handled physically within a station they are 
handled in teletype tape or chad tape form. Hence in 
further discussion when we speak of messages being handled 
in a station it is implicitly meant that such messages are 
in tape form unless otherwise stated. 
1. World Wide Teletype Communications 

The Defense Department World Wide telecommunications 
system controlled by the Defense Communication Agency pro- 
vides for long-haul and point-to-point telecommunications 
between communications stations of the Armed Forces which 
include the NAVCOMMSTAs and other stations within the 
NAVCOMMAREAS. This system is known as the AUTODIN network 
which is a digital high speed message transmission system. 

2. Local Teletype Communications 

The primary means of communications between a 
NENCOMMSTA and» the ships at sea is a teletype system which 
Operates at the rate of 100 words-per-minute (wpm), where 
one teletype word is equivalent to 5 teletype characters. 

There are essentially three methods of teletype 
communication between the ships at sea and the shore based 
NAVCOMMSTA. 

The first is a ship-shore telecommunication link 
which are time shared circuits used by several ships at sea, 
these αἱ τοι ες provide the capability to send messages to a 


NAVCOMMSTA for further relay and delivery. 
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The second method is also a ship-shore link which is 
a dedicated circuit, sometimes called a full period termina- 
tion circuit, between one ship and the NAVCOMMSTA. This 
type of circuit is normally used by the larger ships such 
as aircraft carriers and cruisers, which normally have 
staffs embarked. The full period termination circuit pro- 
vides the capability for these ships to send messages to the 
NAVCOMMSTA for further relay and also to receive from the 
NAVCOMMSTA messages which are addressed to that particular 
ship or embarked staff. 

The third method is the fleet broadcast system. 

This system is the primary means of delivering traffic to 

the fleet as a whole because it provides for simultaneous 
delivery of messages to the operating forces dispersed over 

a large geographical area. This method also employs a 

100 wpm teletype system which could be described as a one 

way system because the broadcast system enables the NAVCOMMSTA 
only to send messages while the ships can only receive 
messages. 

The fleet broadcast system is composed of several 
different types of broadcasts. These are the fleet multi- 
channel broadcast, on which this study is based, a single 
channel broadcast for ships that are not equipped to copy 
the multi-channel broadcast, the submarine broadcast and a 
merchant ship broadcast. Іп some NAVCOMMAREAs these broad- 
casts are all operated by one NAVCOMMSTA or Broadcast Control 


Station (BCS) while in other NAVCOMMAREAs these separate 
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broadcasts are operated by two or more NAVCOMMSTAs. Also 
if the capability of flexibility exists, a separate addi- 
tional broadcast may be created within a NAVCOMMAREA to 
handle special traffic load situation. Sometimes a NAVCAMS 
may shift a braodcast to a NAVCALS in an emergency or to 
Maintain proficiency among the capable NAVCALS in operating 


the broadcast. 


C. MESSAGE PRECEDENCES AND SPEED OF HANDLING 
Naval ción: messages are assigned one of four 
precedences. Each message precedence has a Speed of Service 
Criteria assigned by reference 9, which is a time standard 
for a message to take in transmission from the origination 
of the message to the receipt of the message by the addressee. 
These criteria are shown in Table XVI below. The letter 


in parenthesis is the designator of the precedence. 


TABLE XVI 


SPEED OF SERVICE CRITERIA 


Average Speed 


Precedence "n Limits 
of Seriice 
Flash (2) As fast as Same 
humanly possible 
Immediate (0) 30 minutes 30 minutes-l hour 
Priority (P) 3 hours 1-6 hours 
Routine (R) 6 hours 3 hours-to the 


start of the next 
business day 
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ΤΕ ου τε MULII-CHANNEL BROADCAST SYSTEM 


There are normally five fleet multi-channel broadcasts 
operated in the Naval Communications System. Each broadcast 
which serves a particular ocean area is operated by a desig- 
nated NAVCOMMSTA or Broadcast Control Station. These broad- 
casts have designations of NMUL, EMUL, KMUL, FMUL and GMUL, 
where FMUL serves the Eastern and Mid-Pacific Ocean area and 


is normally operated by NAVCOMMSTA San Francisco [10]. 


A. CHANNEL ALIGNMENT AND USE 

The fleet multi-channel broadcast is a multiplex system 
which has 8 separate 100 wpm teletype traffic channels 
which are transmitted on a single carrier frequency. The 
alignment of the traffic channels is such that the overall 
broadcast traffic load to the fleet is segregated on the 
basis of ship type where each subscriber is served by two 
traffic channels, a primary channel and a secondary channel. 
An example of the overall broadcast channel alignment, 
Channel designations and uses of each channel is shown in 
Table XVII below. 

Only channel pairs 1 and 2, 3 and 4, and 5 and 7 will SS 


considered for this study. 


КО THE RECEPTION SIDE OF THE BROADCAST 
The fleet multi-channel broadcast system for each channel 
is a guard system. This means that each subscriber copying a 


particular channel screens each sequentially serial-numbered 
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TABLE XVII; 


FMUL-BROADCAST CHANNEL ALIGNMENT 


Channel Designation Use 
1 FASW Destroyer force primary traffic channel 
2 FSPC 1 hour delayed rebroadcast'of channel 1 


or overload as necessary 


3 FALD Service, Amphibious and Mine force 
primary traffic channel 


4 FUSN 1 hour delayed rebroadcast of channel 3 
or overload as necessary 


5 FNSC Major warships (carriers, cruisers, 
large destroyers, including flagships) 
primary traffic channel 


6 FOPI Special traffic channel (Not of study 
interest) 

7 КЕЛЕС” l hour delayed rebroadcast of channel 5 
or overload as necessary 

8 FMET Weather traffic channel (Not of study 
interest) 


message sent on that channel to determine whether that 
particular message is addressed to him or nct. If the mes- 
sage is for the ship or for the embarked staff the message is 
processed. If it is not for that particular ship or embarked 
staff the message is ignored and filed. Hence a subscriber 
must perform two necessary actions.’ He must determine from 
the message heading (which contains the addressees) whether 


the ship or embarked staff is an addressee or not. If the 
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message is addressed to the ship then it must be copied with- 
out errors from beginning to end. That is the addressee 

must have a clear enough copy that he is able to obtain all 
the information transmitted without any doubt as to what 

the message says. 

If a subscriber has missed part or all of the heading 
such that he cannot determine whether the message was 
addressed to him or not, or if he has missed part or all of 
the text then we will refer to such a message as a "missed 
message." 

A message may be missed for various reasons, some of 
which are: environmental or atmospheric conditions, which 
cause deterioration of the radio signal, equipment failure 
on the ship or in some Cases error in judgment of the 


operator monitoring the teletype printer. 


C. PROCEDURE FOR OBTAINING MISSED MESSAGES 

If a subscriber copying a particular channel has missed 
one or more messages he must determine whether the message(s) 
was addressed to him and if necessary he must obtain a 
complete copy of the message. He is first encouraged to 
obtain a copy from another ship copying the same broadcast 
channel, providing the subscriber has a separate means of 
communication with the other ship. When the subscriber 
exhausts all local means of obtaining this information his 
next action is to originate a message which is called a 
"service or screen request message." Such a message is 


usually sent to a NAVCOMMSTA requesting them to verify or 


109 





screen their copy of the broadcast to determine whether the 
message(s) in question is for the subscriber. The NAVCOMMSTA 
will in turn reply by indicating which messages are of no 
concern to the subscriber and by rerunning on the primary 
channel those messages which were addressed to the sub- 
scriber. A subscriber which misses only a portion of a 
message may service only for that portion of the message to 
which the NAVCOMMSTA must in turn originate a service 
message in reply giving the subscriber that portion of the 
message missed. 

It is noted that in some NAVCOMMAREAs, the NAVCOMMSTA 
to which a screen request is sent may not be the Broadcast 
Control Station. In these cases designated NAVCOMMSTAs in 
the NAVCOMMAREA also copy the broadcast and serve as 
stations which handle screen requests for a particular 
channel of the broadcast. When these NAVCOMMSTAs process a 
screen request their replys are relayed via AUTODIN to the 
Broadcast Control Station which in turn sends the reply on 
the appropriate broadcast channel. We note that these replys 


are handled in the same manner as first run messages. 


D. HANDLING OF SERVICE MESSAGES AT A NAVCOMMSTA 

Some NAVCOMMSTAs will handle screen requests for more 
than one broadcast channel or in some Cases for more than 
one type of broadcast system. At each NAVCOMMSTA all such 
broadcast screen requests are usually handled at one broad- 


cast service desk. The broadcast service desk handles these 
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Pequests—,on the highest precedence first and on a first in 
first out basis within precedences. 

These service messages are ΜΝ referred to as screen 
requests, since the service desk clerk is required to "screen" 
the heading of the missed messages in the printed copy form 
of the broadcast channel log to determine if the message is 
for the requesting subscriber. Additionally, to clarify 
terminology, each message looked up is referred to as a 
"screen." Thus a "screen request" may contain one or more 
"screens." 

After the service desk clerk determines the result of a 
screen request he then must make teletype tape copies of the 
messages to be rerun, if necessary, which are forwarded to 
the appropriate broadcast channel position for transmission. 
At some NAVCOMMSTAS, when the service desk clerk has finished 
the screen and before he prepares the rerun tapes he prepares 
a short message to the requester indicating which messages 
are of no concern to the subscriber and additionally states 
that those messages which are of concern will be rerun later. 
Other NAVCOMMSTAS use the procedure where they include this 
information with the first message rerun and only after the 
entire screen request has been processed. 

It is apparent that this process is a very time consuming 
task and explains the reason why subscribers are encouraged 
to seek es sources for verification before asking the 


NAVCOMMSTA. 


LEL 





E. TRANSMISSION SIDE.OF THE BROADCAST 

Considering now the transmitting side of the broadcast 
we look at the NAVCOMMSTAs internal handling of messages 
which are destined for transmission on one or more of the 
broadcast channels. Messages arrive at the broadcast area 
in teletype tape form and are, depending on individual 
NAVCOMMSTA configurations, primarily delivered to each broad- 
cast channel position by two 850 wpm reperforators called 
"burpees." The traffic delivered by the burpees is traffic 
which has been received by the NAVCOMMSTA from the AUTODIN 
network or from 100 wpm teletype circuits which terminate at 
the NAVCOMMSTA from any number of sources. The burpees are 
fed from these various sources through a small computer 
system called the Multiple Address Processing Unit (MAPU). 
This system reads the header of a message which contains 
address codes as to where a message is destined within the 
NAVCOMMSTA and routes that message to the appropriate relay 
positions by delivering a teletype tape through the burpees. 
In addition to the burpee deliveries some NAVCOMMSTAs hand 
deliver screen request answers, reruns and other messages to 
the broadcast position. The screen request answers and 
muons are then usually put at the head of the line or queue 
for transmission. 

Messages delivered to the broadcast position are put in 
a grid (queue) by precedence and in order of arrival within 
each precedence. The messages are then taken from the grid 


and put in the channel transmitter distributor (TD), a 
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device which reads the teletype tape for transmission. The 
order in which the ESOS are taken from the grid for 
transmission is on the highest precedence waiting and first 
Mmi first out (FIFO) within а precedence basis. This process 
at some NAVCOMMSTAs is done by one man physically performing 
the queueing process while another man pulls the tapes from 
the grid and places them into the TDs for transmission. At 
other NAVCOMMSTAs, properly equipped, this entire process is 
handled by a small computer system, known as the "tight 
tape" system. In either case the queueing and transmission 


flow is the same for the purpose of this study. 
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APPENDIX C 


GOODNESS OF FIT DATA ANALYSIS COMPUTER PROGRAM 


T. IPES RIÈTION OF CAPABILITY 


The Goodness of Fit Data Analysis Program accepts any 
number of data samples of varying sizes and performs 
prescribed calculations and produces various graphs and 
tables which are used to analyze the fitting of an empirical 
distribution with a selected theoretical distribution. The 
program is written to fit the empirical data to two different 
distributions. These two distributions are the Exponential 
and the Hyperexponential probability distributions in either 
their regular form or the right shifted form. 

The main program uses five different subroutines which 
are named SORT, STAT, HISTO, PLOTP and UTPLOT, the last two 
subroutines being used in conjunction with each other. The 
main program with the subroutines first sorts the data in 
ascending order and prints out the data in this form. Then 
the program calculates statistical parameters and prints out 
the parameter values which are later used in subsequent 
calculations. The next feature provides for sorting the 
data into frequency classes and prints a histogram of the 
data according to the number and widths of intervals as 
prescribed by the programmer. The program then calculates 


ШІ plots the log of the tail distribution for both the 
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empirical and theoretical distributions being tested. A 
description of these calculations along with other mathe- 
matical descriptions are more fully explained below. The 
respective Cumulative Distribution Function (CDF) values are 
calculated and the absolute difference of the distribution 
points is determined and printed. The absolute difference 
provides the user this information for use in performing 
the Kolmogorov-Smirnov (K-S) Goodness of Fit Test on the 
data. The cumulative distribution functions are plotted to 
facilitate a graphical interpretation of the fit of the 
distributions. Finally the data is analyzed by performing 
the calculation of the Chi-Sgquare Statistic and printing out 
a table of results for use in performing the Chi-Square 
Goodness of Fit Test. | 

Features of the program provide a flexibility of 
selecting only certain calculations to be performed and/or 
certain outputs to be printed. These features of selection 
are fully documented in the computer program shown in this 
Appendix. 

This Appendix is not meant to be a user's guide for the 
program but rather to explain and reference the mathematical 


calculations used in the analysis program. 


II. THE THEORETICAL DISTRIBUTIONS 


A. THE EXPONENTIAL DISTRIBUTION 
The exponential density function, cumulative distribution 


function and the associated parameters are given below. 
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f(x). = A ¿X>a, where a is the right 
shift parameter. 

= 0 ,х<а, 
F(x) = l- =A (x-a) ‚x>a, 

= (0 Маса. 
H = E(X) = + а, 
е ж 

р-а 
c^ = VAR(X) = — . 
À 


B. THE HYPEREXPONENTIAL DISTRIBUTION 

The hyperexponential distribution shown below is a 
mixture of two exponential distributions. The observations 
originate in two populations only, with probabilities p and 
1-р апа parameters 2Ap and 2A(1-p), respectively [11]. 

The hyperexponential density function, cumulative dis- 
tribution function and the associated parameters are given 
below. | 


2. -2pA(x-a) 


f(x) = 2p Xe + A х>а, 
- 0 , х<а, 
F(x) = Ee (1-p) LL] , x>a, 
0 , х<а, 


where a is the right shift parameter. 


H = E(X) = > + a, 


Da 
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III. STATISTICAL AND GRAPH CALCULATIONS 


A. SAMPLE STATISTICAL PARAMETERS 

The analysis program has a feature in which preselected 
values for the mean, the right shift parameter a and the 
hyperexponential parameter p can be read in instead of having 
the program calculate these estimates from the data. When 
this feature is not used some of these parameters are cal- 
culated by the subroutine STAT using the familiar formulas 


shown below. 


N N 2 
~ = 2 = sm E 
XBAR = X = > Х., VAREST = N-I > (X, NE: 
i=). Iel 
N 
We note that 
OE: „ы |... 
һепсе, 
LAMEST = bm 
(XBAR - A) °` 


The capitalized names are those used in the program for 


these various parameters. 
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The estimate of the hyperexponential parameter p is 
calculated using the same formula shown above except that 
XBAR = u and VAREST =o. 

The right shift parameter a can also be specified to be 


the minimum order statistic of the data sample. 


B. A GRAPHICAL PROCEDURE OF TESTING FOR EXPONENTIALITY 

When comparing empirical data to the exponential distri- 
bution, for goodness of fit, this graphical technique 
provides an appealing method of either supporting or denying 
the more purely statistical comparisons to be described 
below. 

If the empirical data is really απ ο ρωσ exponentially 
with a cumulative distribution function as shown above then 
the following graphical comparison can be made [12]. For 
the exponential distribution the log of the reciprocal of 
the tail distribution (called the log of the tail for 
brevity) is: 


(х-а) 
u 





y = log Du = , assuming u > 0. 


Thus if we plot y against x, we get a straight line with a 
slope of m and the line will cut the x-axis at the point x=a. 
If we have n data points in the empirical data, we can 

sort the n points in ascending order forming the order 
GE CiStics, 1.е., xj < x, < BE S Xa. The empirical 


cumulative distribution function is defined as 


F(x;) E 


B |=- 
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If we let y = log τεστ hen by pletting уы against X; 
the plotted points can be compared to the fitted straight 
line as a test of the exponential assumption. 

It is noted that the program has the option of using 
a preselected mean or using the sample mean estimate in 
calculating the straight line. Thus, this graphical com- 
parison can be used to make an independent estimate of the 
mean by determining the slope of the empirical data. Then 
the data can be rerun with a preselected mean to attempt 
a better fit. 

When the log of the tail distribution is used with the 
hyperexponential comparison we get a decreasing failure rate 
function for the hyperexponential distribution. That is we 
get a concave curve. In this case the graphical comparison 
can also be used to determine how close the two compare. 

This method of distribution fitting can be used when 
empirical data is analyzed for the first time. Through the 
use of the log of the tail graph we can determine how the 
distribution behaves and find out whether it has a decreas- 
ing or increasing failure rate and thus plan further analysis 


on the data. 


C. THE KOLMOGOROV-SMIRNOV GOODNESS OF FIT TEST 

The Kolmogorov-Smirnov Goodness of Fit Test is a familiar 
test and therefore will not be described here [13].. 

It is noted that the program only calculates the absolute 
difference between the F(x) values of the two distributions 


and prints the results. Therefore the user must scan the 
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points to find the maximum difference and then enter the 
appropriate Kolmogorov-Smirnov table to determine the 
results of this test. | 

The Kolmogorov-Smirnov test is regarded by most 
references as a more powerful test than the Chi-Square 
Goodness of Fit Test but it has one drawback, that being 
that the Kolmogorov-Smirnov statistic becomes inflated when 


the data has large groups of points which are of equal value. 


МЕ THE CHI-SQUARE GOODNESS OF FIT TEST 

The Chi-Square Goodness of Fit Test is also a familiar 
test and will be described only briefly [12,14,15]. 

It is well known that the chi-square goodness of fit 
test has several drawbacks. Among them are its large sample 
character and the dependence upon the careful choice of the 
number and the widths of the intervals. Therefore this test 
should be used with caution. For further discussion Gu 
these characteristics see re ence 14. 

It is also noted that the program has the capability of 
either reading in a specified mean or estimating it. This 
will affect the number of degrees of freedom used when 
testing the chi-square statistic against the appropriate 
table values. The chi-square goodness of fit test with an 
estimated mean uses (k-2) degrees of freedom or (k-1) degrees 
of freedom when the mean is preselected, where k is the 
number of intervals (or classes) used in the test, and where 


the x-axis is divided into the intervals xj < x, <..... «χι η’ 
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The chi-square statistic was calculated in the following 


manner. Let о. be the observed number of data points in the 


sh interval and let е, Бе the expected number of data points 


ше ісі 


interval, where e, = np; and where n is the 
number of data points in the sample being tested and p; 15 
calculated as indicated below for the appropriate distribu- 


tion. 


For the exponential distribution: 
x. +1 

B _ -\(x.-a) -A(x.,.-a) 

қ г 35 À (t a) dt ae 1 E Er 


i π.δ. κ-2, 


and 


> _ DA --а) 
p. = / Ле ВЕ a) at =e ὅπ απ οκ... 
"k-1 


For the hyperexponential distribution: 


Xi+l 
p; -| 2p e 2P (e-2, 2 (1p) e72 17923781 ae 
xs 
-2pÀ (x, -a) РА буута) 
= ple -e 
-2(1-р)А(х.-а) -2(1-p)A(x.,.-a) 
+ 1-p)|e 1 Же JI | 


Sl. K=2; 
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The value e, was restricted to the condition of having 
e, > 5 for each interval [14]. If this condition is not 
satisfied then the observations of that interval are added 
to the previous interval and the appropriate integral is 
recalculated. This is repeated until this condition is 
satisfied. 

Finally the chi-square statistic is determined by this 


formula 


k 2 
X= < ف‎ 
e. 
stat = 1 


This statistic is then compared to the appropriate 


chi-square table value to determine the results of the test. 
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APPENDIX D 
TABEES OF CHANNEL SERVICE STATISTICS FOR MAY-JUNF 
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SIMULATION MODEL COMPUTER PROGRAM 
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